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The design of a centrifugal impeller must yield blades that are aero- 
dynamically efficient, easy to manufacture, and mechanically sound. 
The blade design method described here satisfies the first t wo criteria 
and with a judicious choice of certain variables will also satisfy stre.ss 
considerations. The blade shape is generated by specifying surface 
velocity distributions and consists of straight-line elements that con- 
nect points at hub and shroud. The method may be used to design 
radially elemented and backward-swept blades. The background, a 
brief account of the theory, and a sample design are described in this 
paper. 


The design of impeller blades is a difficult task due to the fact that, 
ideally, three separate and equally important requirements must be 
satisfied. The first requirement is that the impelltT provide acceptable 
distributions of relative velocity on both driving and trailing surfaces of 
the blades in order to minimize the possibility of flow separation and the 
accompanying loss in performance. In addition, the selected blade shape 
must be such that it can be manufactun'd accurately and economically 
by moans of automated fabrication procedures. Finally, the blades should 
be designed so as to keep the stresses at a safe level, eliminating the 
possibility of excessive distortion or fracture occurring during operation. 

For many applications, these three requirements are somewhat in- 
compatible. For example, it is possible to define a blade shape which 
theoretically would give ideal velocit}' distributions throughout. How- 
ever, it is often the case that such a shape is virtually impossible to fabri- 
cate in any reasonable manner. Therefore, any practical design method 
must involve some reasonable compromises with regard to the conflicting 
requirements. The method described in this paper, the prescribed loading 
method, is intended specifically to produce impeller designs which satisfy 
the aerodynamic and manufacturing requirements mentioned above. 
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Although the mi'thod dni'S not include .stress considerations explicitly, an 
intelligent choice of in]nits to tin- method and a careful examination of the 
results obtained will enable' safe de.signs to be produced with a moderate 
amount of effort. 


AERODYN AMIC CON SI DER ATION S 

Lossi's in a centrifugal comjin's.sor are generally caused by boundary- 
layer and S(‘parated flow developments in th(' impc'lh'r and tlu'ir sub- 
seepient effects in the diffuser. Th(' development of the' boundary layers 
and separated flows is directly a result of the velocity distributions along 
the blade, hub, and shroud surfaces. Tfw're are four imjiortant jiarameters 
associated with th(> velocity distribution that affect the boundary-layer 
behavior. They are 

(1) The streamwise velocity reduction. This velocity gradient is 
also recognizc'd as being tlu' main variable that affects boundary-layer 
behavior in two-dimensional flows. 

(2) The ri'duction and increa.se in turbuh'nce intensity along the' 
suction and prc'.ssun' blade' surface's, re'Sjie'ctively. This e'nhiine'es fleiw 
se'jiaratiem alemg the' suctieni surface due' te> jiartial laminariziition. Thee 
e'hange in turbuh'nce inte'nsitv is direectly re'late'd fee the' eliffe're'nce's in 
suctieen and pre'ssure' surface' ve'leecitie'S or leiaeling (re'f. 1). 

(3) The influx of fh)w along the- blade surface's from eene' stre'amtube' 
(a thin annular channe'l ceentaining the; jirimary fleew) tei the' ne'xt cause'd 
by se'condary fleiw. This fleiw usuidly cem.sists e»f leew-meeme'utum fluid; 
the' beeundary laye'rs are' meere' su.sceptible te> se'paratiein whe'u the'y contain 
a large amount eef leiw-meeme'ntum fleiw. The difi’e're'iitial leiading e>f enu' 
stre'amtube' ceunpare'd te> the' ne'xt is dire'ctly rclate'd te> the devcleipment 
e)f secondary fleiws (re'f. 2). 

(4) The highe'r fleew leesse's of the shreeud stre'amtube caused by 
clearance' e'ffe'cts. The'se' le>.sse's cau.se> a pre'.ssure dreep which leads te> fleiw 
separatiem. Incre'ased weerk ininit will offse't the rise' in pressure hisses 
(rc'f.3). 

hi.xtensive expe'rime'iital re'se'arch is curre'utly be'ing carried em to measure 
these effe'cts and re'fine reh'vant theeirie'S. 

The' method .shemld the're'feire' be' cajiable of ge'iie'rating a blade surface 
feir which the following ite'ms are spe'e'ifie'd: 

(1) The streamwise vehicity di.stributiem. This quantity is contreilled 
by both the hub and shremd ceintours and the' blade loading. 

(2) The differe'iice in suctiein and pressure surface ve'hicity (eir load- 
ing) to minimize the turbuh'nce attenuation at the suctiein surface. 
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(3) The difference between the work addition to a streamtube near 
the hub and that near the shroud to minimize the transport of low- 
momentum flow by the secondary flow phenonK'Hon. 

(4) The higher work input into the flow in the shroud streamtube to 
offset the clearance losses. 

From an aerodynamic point of view, it is thus desirable to specify a 
velocity distribution along the pressure and suction surfaces of the blade 
such that the generation of .separated flows is kejit to a minimum. This 
gc'iieral principle would lead to a blade defined by a set of streamlines 
whose individual loading schedule is controlled by the designer. 


MANUFACTURING AND GEOMETRICAL CONSIDERATIONS 


The ideal blade which would be composed of an arbitrary number of 
individually loaded streamlines repri'si'iits a very gimeral three-dimen- 
sional blade surface. It is questionabh' wlu'ther such a blade could be 
economically produced with the nviuin'd accuracy. Impeller casting 
patterns are very difficult to manufacture, particularly when the surface 
of the blade has an arbitrary shapi', when close tolerances have to be 
maintaiiu'd, or when the impi'ller size is small. A more de.sirable blade 
shape is om; defined by .straight-line elements; such a shape; facilitates 
setting up the pattern and allows close- teele-rance-s te> be maintaine-d. 

The machining of an arbitrary blade surface would generally be im- 
])ossibIe be-cause e>f the extre-me- loced curvature of the blade surface. At 
best, the- tool would have peiint ceaitact endy and teee) much time weiuld be; 
reejuire-d tei machine away all the fille-r material. A more aiipropriatc 
machining metheid is eene wlie-re the cutte-r has line ceentact and hence large 
cutting rate. However, the blade shape must still be; ri-stricted to moderate 
inclinations against normals to the hub surface; to avoid tool interferences. 

In summary, from manufacturing considerations, a blade shape con- 
sisting of straight-line elements on both surfaces is preff-rred. The thick- 
ness distribution of the blades is then neces.sarily restricted to a constant 
taper between hub and shroud, and the camber line surface of the blades 
will also consist of straight-line elements. 

Since the camber-line surface is defined by a set of streamlines with 
controlh'd individual loading, the number of streamlines for which loading 
can be controlled and still yield a straight-line element camber surface is 
limitc'd. In general, the individual streamlines form curves in space and a 
straight line can only connect a maximum of three curves in space. Con- 
seciuently, a straight-line element blade can be constructed when the 
loading of three streamlines is specified. These streamlines are assumed to 
be those adjacent to the shroud and the hub and an intermediate one. 
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This specification is considered adequate aerodynamicall}', and the ruled 
blade surface is accepted as a balanced compromise. 

There are three general categories of blades which can be defined by 
straight-line elements. These are termed radial, axial, and arbitrary 
blades. The first two blade types are characterized by the fact that the 
angular orientations of the line elements defining their camber surfaces 
are fixed by the nature of the blades; that is, in the radial direction and 
axial direction, respectively. As a result, the angular orientation (6) can 
be prescribed along only one streamline; the obvious choice is the shroud 
streamline where tlie flow conditions are critical. 

The arbitrary type of blade is characterized by the general and usually 
unspecified orientations of the line elements (excc'pt those at the leading 
and trailing edges which are usually radial elements) defining the camber 
surface. The procedure followed in designing an arbitrary blade involves, 
in general, the loading of the three streamlines in an appropriate manner. 
The orientation of each blade element is then determined by requiring 
that the straight line pass through each of the three space curves defined 
by the loaded streamlines. 

In some cases when for an arbitrary blade the loading is specified at 
three streamlines, the resulting blades may still exhibit a large variation 
in circumferential angle between the hub and shroud streamlines which 
would be unsatisfactory from manufacturing or stress considerations; 
that is, the blade would have excessive lean throughout its length. This 
difficulty was recognized during the development of the procedure and 
special provi.sions were introduced in the procedure to circumvent this 
by specifying at the outset the desired distribution of line element orienta- 
tion and by prescribing the loading along the shroud streamline onh'. 


REQUIRED AND GENERATED INFORMATION 

The prescribed loading method can be utilized to design an impeller 
when its rotative sjK'ed, flow, avi'rage efficiency, overall geometry (that 
is, hub and shroud coordinates, distribution of blade thickness along hub 
and shroud, location of leading and trailing ('dgc, number of blades, and 
exit blade angles) , desired swirl distribution at inlet, and angular momen- 
tum distribution at discharge arc given. The designer must also specify 
any one of four possible types of loading distribution: 

(1) Linear loading from inlet to discharge (fig. la). 

(2) Double linear loading from inlet to discharge (fig. 16). 

(3) Double linear loading from inlet to discharge, maximum loading 
near leading or trailing edge (fig. Ic) . 

(4) General loading distribution (fig. Id). 
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a— Linear loading. h.— Double linear loading. 




('. — Inverse double linear loading. d . — General loading. 

Figure 1. — Available options of loading distributions. 


The choice of a preferable loading is difficult to make, since practically 
no systematic experiments have been performed to determine the relative 
merits of possible loadings. It is suggested that in mo.st applications either 
a constant linear loading distribution, thereby minimizing secondary 
flows (ref. 4), or a double linear loading, minimizing the streamwise 
velocity reduction, will be quite satisfactorj-. Figure 2 shows these two 
loadings superposed. 

The results of the calculation are a full description of the blade surface 
in cylindrical coordinates and the associated thickness distribution. 


LIMITATION OF THE BLADE DESIGN METHOD 


The prescribed loading method can be used to design radial, axial, and 
“arbitrary” blade shapes whose surfaces are defined by straight-line 
elements. However, there are a number of limitations incorporated in 
this approach. These limitations arc discussed briefly below. In some cases 
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Fkhire 2 . — Linear and double linear 
loading superposed. 



tlicse limitations arc' not serious; in others some additional effort is 
necessary before the generated blading can be produced. 

Thickness Distribution 

The gc'uerated blade sliajies an' restricted to having trapi'zoidal thick- 
ness distributions; i.c'., con.stant taper along any liiu' element between 
hub and shroud. An exponential thiekni'ss distribution may lie preferred 
when blade strc'sses arc an imjiortant design considc'ration. This alteration 
may bo made by adding or subtracting thickne.ss from the generated 
surface; no appreciable change in ac'rodynamic variables is expected when 
thickness changes are kept moderate. 

Blade Shape Restrictions 

One fairly important limitation of the pn'.scribi'd loading method as 
applied to arbitrary blades is that it is possible' to obtain unusual shape's 
which may be' difficult to manufacture or which weeuld not be' e'xpecte'd to 
be capable of withstanding the operating stre'sses. Sucli a blade shape may 
result from different heading requirements at the hub and .shroud stream- 
tube'. This blade' can be re'designe'd by spe'cifying a new and more desirable 
straight-line eerientatieen. Feer this case' the' blade heading at the shroud is 
still specified as be'feere, but now the blade' loading at the' hub and the 
interme'diate stre'amline can no longer be' .spe'e'ifie'd. The sample case is an 
illustration of this jerocedure. 

Leading Edge Geometry 

The method gene'rate's the blade .shape' near the h'ading e'dge such 
that the flow ente'rs at zeuee ineide'uee'. Heeweve'r, whe-n the' re'lative inlet 
Mach number is high, this may neet re'sult in an eejitimum blade shape. 
Therefore', se'parate' ])roee'dures are necessary to de'sign a blade h'ading 
e'dge capable of efficiently handling supe'r.soiiic flows. 

Hub and Shroud Contours 

A comph'te aerodynamic impelh'r passage de'sign me'thod would include 
a provision for automatically adjusting the hub and shreiud contours to 
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give a linear or any other prescribed distribution of average relative 
velocity. Due to the phiTiomenon of slip, the average relative velocity 
near the trailing edge decreases and then rises. In order to straighten out 
this distribution, the blade height of the imiieller near the tip becomes 
larger at higher radii, which is usually unacceptable. Considerable effort 
has been expended to incorporate this design feature, but it was found that 
hub and shroud contours which yield a linear relative velocity distribution 
are in general not acceptable, and this hub-shroud contour generation 
procedure was abandoned. 

Trailing Edge Flow 

The spi'cified blade loading does not vanish near the trailing edge as is 
required from the Kutta or “slip” conditions. However, apjiropriate steps 
are taken in the analysis to satisfy the Kutta condition, thereby over- 
riding the .specified blade loading over the last 10 percent of the blade 
chord. In addition, the procedure adjusts the loading such that the 
specified amount of work will still be imparted to the flow. 

GENERAL THEROY OF PRESCRIBED LOADING METHOD 

The overall procedure used in designing an impeller by the method of 
prescribed loading involves the following four basic steps: 

(1) An initial estimate is made of the meridional velocities in the 
impeller passage as describi'd below. 

(2) The meridional velocities are then used in loading the required 
number of streamlines in the prescribed manner. 

(3) A blade surface the shape of which conforms to loaded stream- 
lines is developed. 

(4) A detailed analysis of the flow field in the impeller with the 
devi'loped blade surface is performed. 

Th(> procedure, starting with th(> loading of the streamlines, is repi^ated 
with the new estimates of the flow field obtained in the detailed analysis. 
Experience has shown that two repetitions of the procedure are sufficient 
to ensure convergence. This is defined as occurring when the estimated 
flow fli'ld used in loading the streamlines and defining the blade shape is 
approximately the same as that obtained in the .subsequent detailed 
analysis. The various steps involved in the design procedure are discussed 
below in more detail. 

INITIAL VELOCITY ESTIMATE 

The di.stributions of meridional velocity along the shroud, hub, and 
intermediate streamlines arc required for use in the loading procedure. 
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Initial estimates of these distributions are obtained by means of the one- 
dimensional analysis. First, a number of normals are constructed which 
are normal to both the hub and shroud at the points of intersection. These 
normals to the hub and shroud are circular arcs. 

Neglecting curvature effects, the mean meridional velocity at a normal 
is given by the continuitj’ equation 


C 

^ m,m 


m 


( 1 ) 


The density p„ is a function of the mean meridional velocity, C^.m, and 
equation (1) is solved by an iterative proce>dure u.sing ce)inpre.s.sible; flenv 
equations to relate static density to total temperatures anel pressures. 
The values of C„,„, calculated in this manner are taken as the initial 
estimate's eif the meridional veleicities along the interme'diato streamline. 

If it is a.ssumed that zerei circulation exists in an are-a beninded by two 
adjacent normals and the hub and shroud streamlines, the following 
expression relating (\,, and can be derived. 

Cm.h dmh = Cn,.s dnis (2) 

If the additional assumption is made that the quantity pmCmf varies 
linearly with normal le-ngth be'tween hub and shroud, the'ii feir ceenstant 
density 

^ m ,mPm 2 m m (H) 

A rearrangeme'iit and simultaneeous seilution eef e'quatieins (2) and (3) 
results in the folleiwing e'xplicit relations for r„,„, and 


(r,/ri,) + {dms/dmh) 

(4) 

r 

^ m ,h , ^ m , s 

dnih 

(5) 


The derivatives dm^/dnii, are obtained from thc' ratio of mi'ridional length 
between two normals at hub and shroud. The initial e.stimatc's of the 
meridional velocities at the calculation stations along the hub and shroud 
are then obtained from the latter two eipiations along with the inter- 
mediate values, Cm,m- 


STREAMLINE LOADING PROCEDLRE 

At the outset of the proct'dure for loading a streamline, the mi'ridional 
coordinates (r,z), figure 3, of a number of points or calculation stations 
along the streamline and the values of meridional velocity (Cm) at these 
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MERIDIONAL PLANE 

P’loiiKE 3. — Xomenclaiure of relative flow and hladc-elemeni geometry in a centrifugal 

compressor. 


jioilits are known from thn n'sults of the previous steps in tlie design 
method. First, the distribution of angular momentum is calculated from 
the given loading distribution IF,— 11’^ from 


d{rCe) 

ds 


Z 


(Tr,-ir,) 


(6) 


where s is the camber-line length and Z the number of blades. 

Assuming for the moment that the blade camber-line distance from 
the leading edge is also known at each point, then the swirl can be dis- 
tributed in the specified manner and a value of rCe can be computed anj’- 
wherc along the streamline'. Combiiu'd with the meridional velocity dis- 
tribution, sufficient information is then available to construct the velocity 
triangle and evaluate the flow angle (/3) at each calculation station from 

/ Ce~rui\ . 

(“ct) 

In the above description of the loading procedure, it was assumed that 
the blade camber-line distance from th(> leading edge was known at each 
calculation station. In actual fact, this is not the ca.se at the outset of the 
proci'dure. Therefon', this (juantity mu.st be determined iteralively. The 
initial estimate of the blade camb('r-line distance is taken as the meridional 
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arc length from the loading edge. The streamline is then loaded in the 
manner described to obtain values of /3 along the streamline. The angular 
coordinate (6) of the blade at each station is obtained b\' integrating the 
calculated /^-distribution as follows; 

/•'"tan/3 , 

6= / dm (8) 

T7i 1 ^ 

In this iiit('gration, 0 is arbitrarily initializc'd to be zero at the trailing 
edge point on tlu' shroud. Finally, the three-dimensional coordinates 
(r,z,6) of the blade cambc'r line are u.s(>d to eomi)ut(' the camber-line 
distance at each point. Th(> procedure is rc'peated using th<' calculated 
values of camber-line distance as the new estimates. This is done until 
tlie estimated and calculated values of 1h(‘ total blade camb('r-line distance 
from leading edge to trailing edge are tlu> same. An additional iteration 
on the distribution of angular momentum is neces.sarj- when geometrical 
constraints require the relative streamline to pass through a given angular 
interval. Such conditions arise whenever the blade surface is to be deter- 
mined from th(' loading specification on more than oiu' streamline. The 
requirement follows from the stipulation that all streamlines int(>rsect 
with an initial and final straight-line blade (‘lernent. 

The relative strc'amlinc' coordinates that r(>sult from a specific'd loading 
distribution and the meridional velocity distribution defiiH' the blade 
cambc'r line for an infinite number of blades. A modification of this camber 
line is introduced to account for flow deviation associat('d with a finite 
number of impellcT blades (i.e., slip). For the major ]iart of the blade 
passage, the average flow direction is a.ssumed to coincide with that of the 
blade camber line. In th(‘ latter third of the blade i)assage, a progressive 
degree of deviation is assumed to occur. At the impeller discharge the 
amount of deviation is efpiivalent to the slip. 


BLADE COA STKUCTIOA 

The third major st(‘p in tlu' procedure is the synthesis of a mean blade 
surface from the thr(‘e camber lines that hav(‘ been gcaierated iii the 
pr(‘vious st(‘p and the two line elements at the leading and trailing (>dge. 

The g('ueral procedure' for constructing a line' eh'ment is illu.strated in 
figure 4 which shows thre'c gem'ral space' curves and an initial element 
intersecting each e>f them. The' incre'me'iital arc h'ngtli. As,-, aleeng the first 
curve freem the initial line ('lernent is cheise'n arbitrarily to provide the 
origin of the adjacent e'lement v hich is to be ceenstructed. The incremental 
arc le'iigths, As; and A.s'a-, rejirese'iit the distance's from the' initial eh'ment 
to the points of interse'ction of the adjace'iit e'le'ment with the re'maining 
two curve's. It can be shown that the' distance' ratios de'fining the adjacent 
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line element are given by the following expressions: 



These relations are based on the property of ruled surfaces by which the 
infinitesimal displacement of a straight-line element is related to vectors 
which are tangent to the surface at different points on the line element. 

In the above equations, the unit vector, 1 , represents the direction of 
the initial line element from the first curve, and t is the unit tangential 
vector of each curve at the point of intersection as shown in figure 4. The 
value of t is found from differentiation of the previously defined camber 
line with respect to the distance s measured along it. 

The ])roc('dure for defining the camber surface* of the blade: is begun by 
ge'uerating em the* loaded shroud streamline be-tween the leading and 
trailing e*dge a numbe>r of eciually spaced pehnts representing the origins 
of the line; elements te> be constructed. Starting from the leading e>dge 
eleme'iit, the points e>f intersection of the next line with the intermediate 
and hub streamlines are determined by computing the incremental arc 
lengths along tlu'.se curves with (‘(piations (9a) and (9b) . The construction 
procedure is performed in a .stepwise fashion with the generation of each 
line element being the basis for the gi'iieration of the next one until the 
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trailing (‘dge is reached. This process is neces.sarili’ an iterative one, .since 
initially the exact values of t and 1 are not known. 

As has been point('d out bc'fore, the method for constructing a blade 
surface from thrc'e streamlines in space ma\' lead to undesirable geometries. 
If this occurs, the ca.se may be rede.signed by specifying the orientation 
of the straight-line ehmients in advance in t(>rms of two angles, «« and a,, 
in addition to the loading distribution along the shroud line. 

The straight line and its normal projection on the meridional plane 
include' an angle of lean or dihedral (ae) as shown in figure 3. The angle 
between the projection and the impeller axis defin<>s the slope (angle a^) 
of the straight-line oleme'nt. The describe'd definition of blade element 
orientation leads naturally to the following oftc'ii used s])ecial blade 
surfaci's: 

(1) The radial ('lenient blade surface is determim'd by a camber line 
through which radial line' elements are drawn to define the surface. The 
elenn'iits an' slojicd v('rtically without h'an (a 9 = 0, a^ = 90°). 

(2) The axial ('lenient blade surface is determim'd by axial liiu' 
elements originating from a camber liiK'. These eh'iiK'nts are slopc'd 
horizontally without lean (as = 0, q!z = 0). 

(3) The arbitrary blade surfac(' is defined by straight-line (‘lenn'iits 

('iiianating with arbitrary ori('iitation from a camber line. TIk' (‘lenients 
may have slope and h'an az9^0). 

The camber surfac(' and the given distribution of normal blade thick- 
iK'SS along th(' hub and shroud contour (h'termiin' the suction and pn'ssure 
surface of th(' blad('. At ('ach point of inters('ction of cambc'r line and 
straight-liiK' eh'ment, the normal thickness is inter])olat('d from the giv('ii 
data and one-half of this value is appli('d in each direction of the local 
normal to the mean blaih' .surface to (h'fine points on the' outer blade 
surfaces. The remainder of thes(' out('r surfaces is again described by 
straight-liiu' ('h'ments joining corresponding hub and shroud points. For 
purpos('S of the following flow analysis, th(' blade geoiiK'try is defim'd by 
a matrix of discrete point c()ordinat('s and tang('ntial blade thickiu'.sses 
on the straight-line ('lements of the mean black' surface. The interpolation 
of blade thickn('.s.s(‘s and subsc'ciuent conversion to tangential thickiK'.ss 
recpiirc's the determination of th(' local normal to the mean blade surface. 
On straight-liiK' eh'mc'nts this normal vc'ctor is a liiu'ar combination of 
the known normal vectors at the hub and .shroud points of the elemc'nt. 
During the flow analysis, ])lanar interpolation of the black' geometry is 
used for any blade' point not ('X]ilicitly locate'd on a S])('cifi('d straight- 
line element. 
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DETAILED FLOW ANALYSIS 

The final sti'p in the blade design procedure is an analysis of the flow 
field that corresponds to the di'signed impeller blade passages. The 
analysis serves to reevaluate the meridional flow field from which the 
blade shape has evolved, by means of prescribed loading distributions on 
one or more streamlines. Three cycle iterations of the blade design are 
usually sufficient to achiiwe convergence on the meridional velocity 
distribution. 

The most conveniimt way of expressing the governing equations is in 
so-called streamline coordinates. These coordinates lie in a meridional 
plane (see fig. 3). A set of fictitious streamlines is defiiK'd as the inter- 
sections of annular stream surfaces and the meridional plane. The stream 
surface itself is obtained by rotating the actual streamline for the infinite- 
number-of-blades solution around the axis. 

Using the distance along streamlines (the meridional distance), m, and 
the distance normal to the streamlines, n, as coordinates, the governing 
equations reduce to the following; 

(1) The normal momentum equation, a differential equation which 
primarily determines the distribution of meridional velocity along 
normals. 

(2) The continuit}’ equation, an integral equation which primarily 
determines the position of the streamlines along normals for a given 
meridional velocity distribution. 

(3) The conservation of angular momentum ec(uation (upstream and 
downstream of the wheel) or the Euler ivork equation (within the wheel), 
which is an algebraic equation applied along streamlines betiveen normals. 

(4) The energy equation, an algebraic relation between the change in 
stagnation enthalpy and the energy in]iut along streamlines. 

(5) The equivalent of the streamwise momentum equation, an 
algebraic relation which primarily determines total pressure loss along 
streamlines. 

These equations are generally expressed in two slightly different forms 
depending upon whether the flow region consideri'd is inside or outside 
the blade passage. In addition, for gases, the equation of state provides 
the necessary relationship between fluid properties. The remaining equa- 
tions, including the many useful reflations derivable from the (>quation of 
state and the di'finitions of stagnation propcrtii's, are well known and 
not presented here. Detailed discussions of this flow analysis method may 
be found in references 3, 5, and G. 

Brieflj-, the solution procedure is an iterative one which involves first 
estimating the positions of the streamlines throughout the flow field, 
solving the normal momentum equation and the continuity equation 
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simultaneously along each normal (using other ('quations to rc'late flow 
properties on each strc'amline to those of the upstream normal), revising 
th(' estimate of streamliiu' jiositions, and reiieating the process until con- 
vergence on the stn'amline jiosition is obtaiiu'd. The solution so obtained 
is based on the assumption that thi' flow is everywhere parallel to the 
blade surfaces and, in (‘ffect, accomplishes the determination of appro- 
priate streamtubc' wiiiths throughout the flow field from which the 
meridional vi'locities are di'termiiK'd by the use of (luasi-one-dimcnsional 
flow relationships. 

SAINIPLE mPELLER DESIGN 

The blade design method is illustrated by a sample case with the 
following design specifications; 


Rotational speed 

27 900 rpm 

Design mass flow 

4.61 Ibm/sec 

Inl(*t total pre*.ssure 

14.7 psia 

lnle*t total tempe*raturc 

528° R 

Discharge* total te*mperature 

o 

00 

O 

Impe*ll(*r isentrojiic e*fficiency 

0.945 

Impeller discharge blade angle* 

1 

O 

Impelle*r inlet hub radius 

1.5 in. 

Impelle*r inlet shroud radius 

3.33 in. 

Impelle-r discharge radius 

5.08 in. 

Im])elle*r discharge width 

0.539 in. 

ImpelleT axdal l(*ngth 

2.75 in. 

Leading edge* e*lement orientatieni 

a. = 90° 

o?e = 0 

Trailing e*dge* e*le*me*nt orientatiem 

a, = 0 

«,= -30° 

Blaele number 

16 

Blade* thickne*ss alemg the shroud 

0.05 in. 

Ave*rage blade* thickne'.ss alemg the* hub 

0.15 in. 


The impeller hub and shroud contours arc' shown in figure 5. An ar- 
bitrary mean bladi' surface of unspecified bladi' ek'mc'nt orientation was 
selected for the initial di'sign. Th(' loading distribution along the shroud 
streamline was chosen for minimum diffusion on the blade* suction and 
pressure surfaces. This loading schedule is as .shown in figure* 15. Alemg 
the hub and an inte*rme*diate* stre*amline*, a line*ar loading sche*dulc was 
maintained as shown in figure* bi. At the* blade* in!e*t this loading was ad- 
just e*d to meet the* ge*eim(*trical coirstraint that all thre*e* .streamlines should 
pass through the* initial .straight-line* ele*me*nt at the* le*ading e*dge* and the 
final line element at the* trailing e*dge*. 
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Figure 5. Impeller contours for the datian stage {25-degree hackslope impeller). 


I'igure 6 shows the blade element orientation resulting from the initial 
design. The leading edge blade element is radial while a negative or 
forward lean, as, is indicated in the inducer region. The blade lean reverses 
its sign toward the middle of the blade and has a negative or forward lean 
of 30° at the trailing edge as prescribed. Blade element slope, Uz, varies 
uniformly from vertical at the leading edge to horizontal at the trailing 
edge. 

Several design considerations led to a revision of the blade element 
orientation. A first concern arose from blade root bending stresses at the 
inducer inlet. The subject impeller forms a basic design standard in a 
series of research compressors that represent systematic variations of 
design parameters. For this reason, a simpler variation of blade element 
orientation was selected for the impeller by specifying zero blade lean 
throughout and a variation of slope as shown in figure 6. This design 
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Normalized Camber-Line Distance 


FuiriiK f). — lilailc-eleiimil nrienhilion for xaniple impeller design. 

modification (‘liminatcd blade root bendinp; stri'sscs in the indueca- and 
the originally calculated negative blade lean, ae, at the impeller exit. 

The final velocity di.stributions on the blade surfaces arc' shown for the 
shroud strc'amtube in figure' 7. Thc' velocity variation along the' shroud 
was not affectc'd by the rc'visc'd distribution of blade elc'ment oric'iitation, 
since the loading schedule on the shroud strc'amline was not alterc'd. 
Within tlu' numc'rical accuracy of this apjdieation, the- loading along the 
shroud strc'amtube varies in the' prc'seribc'd manner. A suitable modifica- 
tion of the hub and shroud contour would have' produced a more linear 
variation of the average relativ'c velocity and a rc'semblance to thc' 
idealized velocity di.stributions of figure 1, but the prc'sent distribution 
was considerc!d satisfactory. 

Thc final velocity variation along the hub strc'amtube as shown in 
figure 8 was not significantly affc'ctc'd by the revisc'd orientation of blade 
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Camber-Line Distance, in. 


Fkji’RK 7. — Design-point tilade-surface-velociti/ distritiiilions for soiiiplc impeller {shroud 

strcaiiduhe). 

clemoiits, altliough blad<‘ loading was inproased along tla- forward jiortion 
of tho hub camber line. Th(‘ large variation in average relative velocitj’ 
along the hub could be reduced by modification of the hub contour but 
was considen'd acceptabh' in its jircsent form. 

CONCLUSIONS AND RECOxMMENDATIONS 

The impell(>r blade design method has been utilized to geiu'rate impellers 
with 0°, 25°, 45°, and 55° backslope at the trailing edge. Some of these 
designs were cast and others machined, but due to the straight-line ele- 
ment principle the manufacturing of these designs was relatively straight- 
forward. The test results of these designs showed that their aerodynamic 
performance was very satisfactory. Dynamic and static pressure measure- 
ments have shown that the actual surface velocity distributions agree 
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FioruK 8. — Dexiyn-point bldde-Kiirfare-velocih/ diKtrihulinnn for impeller {huh rtreamlubc). 


closely with those; calculated with the prc'scribed loading ]iroc(‘dure as 
described above. 

The method involves a large and complicated set of calculation stc'ps 
which are not presented here'; the'se ste])s are ceimpute-rized so that a 
solution is readily obtainable. 

There are twee are>as whe're' future activitie's in impelk'r blade design are 
essential. Tlic first is linking the ve'locity distribution along the; blade 
surfaces at different streamlines to the impe'ller and impeller-induce'd 
diffuser lees.ses. When this link is e'stablishe'd, there is more confidence that 
a specified velocity distributie)n will achieve low impeller losses. The 
second are'a is the atte'nuatieui of impe'ller losses by various boundary-layer 
treatments such as suction, bleiwing, split blade's, tandem blades, and wall 
treatment. Both are'as are curre'iitly under investigation in various 
research establishments. 

A logical extension of the prese'iit method is to gene'rate directly from 
the computer program a paper tape to be used in a numerically con- 
trolk'd five-axi.s milling machine for manufacturing impellers which 
cannot be easily cast at present (i.e., titanium) . 
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LIST OF SYMBOLS 

A Angle bctiveen meridional velocity and impeller axis, deg 
C Absolute velocity, ft/ sec 

1 Vector tangent to line element, ft 

m IMass flow rate, Ibm/sc'c 

TO Streamline meridional distance, ft 

H Distance along normal to meridional streamline (positive in 
direction of increasing radius) , ft 
r Radius, ft 

s Camber-line distance, ft 

t Vector tangc'iit to streamline, ft 

t Blade thickness, ft 

W Relative velocity, ft/scc 

Z Number of blades 

z Axial distance, ft 

a Angle of blade surface orientation, d(‘g 

/3 Relative angle from meridional plane (positive in direction of wheel 
rotation); when overscon'd (d), signifies flow angle, deg 
d Cylindrical coordinate angle around s-axis (positive in direction 
of rotor rotation), deg 
p Density, Ibm/fV 

o) Angular velocity, rad/sec 

Subscripts 

1 Impeller inlet 

2 Impeller exit 

d Pressure surface 

h Hub 

m Meridional value 

m,m Mid-passage meridional value 

r Radial direction 

s Shroud 

t Suction surface 

z Axial direction 

6 Tangential direction 
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DISCUSSION 


F. DALLF.XBACH (Consiiltanf , U.S. Army): Th(‘ authors prcsi'iit a 
method whereby tlie l)lade sliapi' may he determined by siieeifving the 
lilade loading distribution along the floM' path. Starting with a prescribt'd 
meridional shape, the m<‘ridional velocities are computed along the 
shroud and hub 1)3- assuming that zero vorticit3' exists along the stream- 
line normals and that the quantit\- pC„,r varic's linearh' with normal 
length between hub and shroud. These' assum]itions are incorrect, since 
the vorticit\- along the normals in general is not zero but varies with the 
normal distance and the normal gradient of the (juantity pC„,r is not 
constant. The resulting meridional vi'locities must then'fore b(‘ considered 
(piite approximate and, in fact, may be in considerable error from the 
correct values. 

The relative flow angle /3 is computed bv prescribing the blade loading 
along a particular streamline as a function of the camber-line length s. 
The angular momentum rCe is then readily calculated. However, until 
the relationship between the camber-line length s and tlie meridional 
length m is established, the angle' cannot be de'termined. 

The solution of the camber-line lengtli as a function of the meridional 
length leads to the solution of a nonlinear differential equation in which 
the overall camber-line length is a variable, d'he final solution which 
ma3' be obtained b3' numerical methods must satisf3’ the given overall 
meridional length. This procedure is verv time consuming. 

If, instead of specif3’ing the blade loading distribution along the camber- 
line length, the angular momentum rCe is specifu'd along the meridional 
length, the angle /3 ma}- be directh' computed and the blade-to-blade 
velocities readil3' determined along the meridional and camber-line lengths. 
This procedure requires little time and affords a readv means of adjusting 
the blade loading in the event the first trial does not give the desired 
result. 

B3- specifying blade loading distributions along the shroud and hub 
contours, the total angular displacement between inlet and exit is de- 
termined from the distribution of the blade angle (not the flow angle). 
As pointed out b\- the authors, the total angular displacement for the 
shroud and hub should be essentiall}' the same in order that the leading 
('dge should be a straight radial line. The iiroblem of specifying blade 
loading distributions for both the shroud and hub contours for backward 
curved impellers is that the angular displacements of the shroud and 
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hub are, in gi'iieral, not (‘([ual. The authors rerognizi' tliis and point out 
that adjustments of the sliroud and liul) blade loadings ar(> required to 
satisfy (‘qual angular disiilaeenu'nt l)ut they do not trivat in any detail 
what comjiromises are lU'cessary. 

The authors’ ])roposal to generate th(> inqiellc'r tilade surfac(‘s by straight 
liiu's connecting the hub and sliroud is cjuiti' jiractical and, iiidei'd, is 
being done today where such impc'llers arc' machined on a five-axis milling 
machine. 

The authors’ discussion of the impc'ller dc'sigu philosophy is very well 
presented and they are to be commended for presenting this paper 
covering a new design approach for cc'iitrifugal impellers. 

D. G. WILSOX (Massachusetts Institute of Technology) : lliis is a 
reciuest for further information about the actual boundary-layer bc'havior 
on the tc'sted impc'llers. It is statc'd that mc'asurc'mc'nts of pressure distri- 
butions have shown good agrc'c'iuc'iit with prc'dictious, and this agrc'cmc'iit 
is very heartc'iiing for iiropcmc'nts of analytical dc'sign mc-thods. It would 
be very heljiful to be shown the mc'asuremc'nts. 

In particular, such a dc'mcmstration of success would help to (luiet a 
strong criticism of thc'oretical methods of dc'sign of cc'iitrifugal comiiressors. 
It is c'asy to show that, unless the radial velocity at outlet is made un- 
acceptably high, thc' mean vc'loc'ity through a cc'iitrifugal compressor 
with nc'ar-radial blades must undc'rgo a decc'leration which in two-dimc'n- 
sional flow would result in sc'iiaration. Conditions must tlierc'fore be worse 
for the pressure surface' than for the mean, and worse' at the hub than at 
the' shroud. Figures 7 and 8 confirm these rathc'r otivious statc'uieuts. 
But flow analj’sis by thc'sc' methods is useless once separation has taken 
place. Do the authors agrc'c? Does separation take place nc'ar (‘iiough to 
the im]ieller outlet that nc'ithc'r the work output (slip factor) nor the 
impc'ller c'fficic'ucy is grc'atly affected? Evc'n so, the prc'.ssurc' distributions 
would be strongly, if locally, affc'cted. Or do the secondary flows, instc'ad 
of causing earlier sc'iiaration as one might imagine, thin down the im- 
portant jiarts of the boundary layer and concc'iitrate all the low-momen- 
tum fluid in one noucritical location? 

I believe that the significance of these quc'stions is that, although many 
previous tests on centrifugal comprc'ssors using flow visualization of 
various types have' shown serious flow brc'akdown which has often been 
ascribed to sc'condary flows, none' of thc'se impc'llers, to my knowlc'dge, 
was designed in thc' first place' to have the kind of controllc'd velocity 
distributions which it has bc'C'ii the aim of the authors to produce. They 
are to be congratulated on a stimulating jiapc'r. 

H. D. LIXHARDT (.Mrco Cryogc'iiics) : The' authors are to be con- 
gratulated for reporting on their continuc-d c'ffort in the dc'veloiunc'nt of 
ac'rodi'namic inqicller dc'sign methods for centrifugal comprc's.sors. Since 
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flip conijilox flow field of ci'iitrifugal and, in ]iarticular, high pressure ratio 
eentrifugal eompri'ssor.s is not completely understood yet and does not 
lend itself to rigorous analytical treatment, approximate solutions and or 
analysis of assunu'd flow fields are being utiliz(>d for most compressor 
designs. Tlu' authors’ suggested loading te'chniepie, which lacks exiieri- 
mental foundation and physical understanding of the impeller flow, falls 
inte) the> latter cate'gory of appreiximate solutions. The ap]ilication to 
practice' eif the authors’ me'theid re'cjuires a large and complicati'd set e>f 
calculatiem ste'ps (the authors neglected to pre-se'iit and discuss the eejua- 
tions) which can only be handle'd by a ceimputer with a large storage 
capacity. I am sure continued experimental and analytical effeirt will 
eventually refine the preeposed technieiue, thus presiding a rea.sonable 
de'sign eif centrifugal compre'ssers e>f mode-rate e>r me'dium pressure ratio. 

Ceensidering industrial applications, erne should ned feerget e>r under- 
e.stimate the classic hydraulic appreiach as enitline'd in re'ference D~l. 
The advantage of the “Conformal Transformat ie)ii Technicpie” is that 
the designer can prescribe the blade leiading acceirding te> his jie'rsonal 
understanding ejf the anticipate'd fleiw and his ex]3e'rience with manu- 
facturing eif similar machine's. In additieen, a cle-ar picture of the impeller 
passage can be readily obtaine'd, which alleews application of basic diffusieen 
and boundary laye'r stability analysis. The ceeiiformal transfeirmatiein e>r 
geometric technique has bee'u applied succe's.sfully in the past to in- 
dustrial pumps and rocket type turbopumps and presents a good and 
rapid de'sign tool feer custom-e'iigineered industrial jiroce'ss compressors. 
As a matter of fact, the conformal tran.sformation tee'hnique has re'cently 
been computerized at Airce) and successfully applied te> medium and high 
pre.ssure ratio industrial ceimpressors. I'or this purpose the computer- 
generated blade coordinates were punched on a jiaper tape which has 
been utilized for automatic machining of impeller jiassages with a tape 
controlled machine. For exanijile, figure D-1 shows a high specific speed, 
mixed flow, 2:1 pre.ssur(‘ ratio, 5000-hor.sepower process compressor that 
achic'ved 86-percent efficiency at design point and exhibited a wide flow 
range duo to the 85° back-swept design. 

.JAXSEX AND KIRSCHXFR (authors) ; Mr. Dallenbach raised three 
points, and thi'y are discussc'd below. 

Th(' meridional velocities calculated by assuming zero vorticity are 
those used as the first values in an iterative scheme. As is ]iointed out in 
the t('xt, subsc'quent velocities in the iteration ari' determined by a de- 
taih'd flow analy.sis of the type described in references 3, 4, 5, or 6. We 
agree that the initial estimates are quite approximate, but the final 
velocities are accurate. 

W(' agree that the procedure for loading the blade along the camber 
line is very time-consuming, since again an itc'rative scheme is rciiuired. 
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FiGriiK l)-l. — r>000-horRepou'er compressor wheel. (Airro Vrijogenics) 

Howi'vcr, the vf'locity distrilnition along the blade length (and not that 
along the meridional haigth) is the item to be eontrolh'd, sinee it govi'rns 
the boundary-layer eharaeteristics and si'condary flows. Thus, a tc'eh- 
idcally superior procedure was preh'rn'd over a method which yieldi'd a 
shorter com])uter running tinu'. 

The compromise which has to be accepted when the blade orientation 
{ae and ad is siiecified is that the velocity distribution along the hub 
and an intermediate streamtube cannot be siiecified. Howi'ver, the shroud 
streamtube velocity distribution should still !)(> specified. 

Professor Wilson asks valid ([uestions about the applicability of theo- 
retical methods to eentrifugal impeller design. What the present paper 
attem])ts is to bridge the gaj) betwc'cn resi'arch and di'sign. From n'si'arch 
w(‘ know that boundary-layer separation mon* than anything else will 
geiu'rate large losses. .\lso from research we have idi'iitific'd the significant 
flow mechanisms that govern separation in im])ellers. Each of these flow 
mechanisms may be controlled by manipulation of tlu> flow field and this 
liapc'r provides a jiroci'diin' by which this may be accompli.shed and 
which yields a producible' impellt'r. It is true' that more ri'search is needed 
to associate om- flow field with a superior im])eller iierformance, but 
current work at various ('stablishments, including the authors’ own firm, 
has uncovered some dc'linite' trends in vi'locity distributions (in addition 
to those associated with limiting the irdet to discharge' re'lative ve'leicity 
ratiei) which improve impeller pe'rfeirmance'. 
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The iisofulnoss of tho proci'duro is further exi'nijilified liy considering 
the correlation between calculatc'd and mi'asun'd results. Figures D-2 
and D-d show the calculated and superposed measured velocities along 
the shroud streamtube for an iniiieller with 50° backslopi' generating a 
jiressuri' ratio of 2.4 at an rpm of 28 000. The im]iellers for figures D-2 
and D-3 are similar in all aspects, except that the exit width for the case 
of figure D-3 is 13 percent less than that for figure' D-2. The' measured 
data shown are reduced from measurements of average static pressures 
and dynamic pre.ssures taken along the shroud, tog('th('r with total pre.s- 
sure traverses at the* impeller discharge'. The agre'cment betwe'en calcu- 
late'd and measured data is rea.sonable' in the' initial 60 jierce'ut eef the' 
blade passage. Afte-r this ]ioint, the me'asured velocitie's are' higher, which 
indicate's that the' drew is separated. The' actual pas.sage flow is compared 
with the calculated flow in figure' D-4. By prescribing a more reaseenable 



FrcirKE D-2. — Coinparii^on of calculalctt ami atfasured hUide-i^urface-velonltj doilnindion 
(dong the i^hroud for oO-degree hacktdope impeller {separation starts at a cainbcr-line 
distance of 2.9 inches). 
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velocity distribution, a new impeller can now be designed with the pro- 
cedure. The occurrence of separation can thus be delayed or eliminated in 
the new impeller, ensuring higher efficiencies. 

The writers appreciate Dr. Linhardt’s comments and interjiret his 
discussion as a brief outline of an alternative blade design method. 
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